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ABSTRACT 

GRS 1915+105 was observed by the CGROIOSSE 9 times in 1995-2000, and 8 of those observations were 
simultaneous with those by RXTE. We present an analysis of all of the OSSE data and of two RXTE-OSSE spectra 
with the lowest and highest X-ray fluxes. The OSSE data show a power-law-like spectrum extending up to > 600 
keV without any break. We interpret this emission as strong evidence for the presence of non-thermal electrons 
in the source. The broad-band spectra cannot be described by either thermal or bulk-motion Comptonization, 
whereas they are well described by Comptonization in hybrid thermal/non-thermal plasmas. 

Subject headings: accretion, accretion disks — binaries: general — black hole physics — radiation mechanisms: 
non-thermal — stars: individual (GRS 1915+105) — X-rays: stars 
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1. INTRODUCTION 

The black-hole binary GRS 1915+105 is highly variable in 
X-rays (Belloni et al. 2000, and references therein). Still, even 
its hardest spectra are relatively soft, consisting of a blackbody- 
like component and a high-energy tail (Vilhu et al. 2001). They 
are softer than those of other black-hole binaries in the hard 
state, which EFe spectra peak at ^ 100 keV (e.g., Cyg X-1, 
Gierliriski et al. 1997), and are similar to their soft state (e.g., 
Cyg X-1, GierUnski et al. 1999, hereafter G99; LMC X-1, LMC 
X-3, Wilms et al. 2001). 

The blackbody component arises, most Ukely, in an 
optically-thick accretion disk. On the other hand, there is no 
consensus at present regarding the origin of the tail. All three 
main models proposed so far involve Comptonization of the 
blackbody photons by high-energy electrons. They differ, how- 
ever, in the distribution (and location) of the electrons, which 
are assumed to be either thermal (Maxwellian), non-thermal 
(close to a power law), or in a free fall onto the black hole. 

A discussion of these models is given in Zdziarski (2000), 
who shows that the thermal and free-fall models of the soft state 
of black hole binaries can be ruled out, mostly by the marked 
absence of a high-energy cutoff around 100 keV in the CGRO 
data (Grove et al. 1998; G99; Tomsick et al. 1999; McConnell 
et al. 2000). The present best soft-state model appears to in- 
volve electron acceleration out of a Maxwellian distribution 
(i.e., a non-thermal process), which leads to a hybrid electron 
distribution consisting of both thermal and non-thermal parts 
(Zdziarski, Lightman & Maciolek-Niedzwiecki 1993; Pouta- 
nen & Coppi 1998; G99; Coppi 1999). 

In this Letter, we present all OSSE observations of GRS 
1915+105. We then choose two OSSE spectra corresponding 
to the lowest and highest X-ray flux and fit them together with 
spectra from simultaneous RXTE pointed observations. The 
spectra, showing extended power laws without any cutoff up 
to at least 600 keV, provide strong evidence for the presence of 
non-thermal Comptonization. More extensive presentation of 
the combined X-ray/OSSE data will be given elsewhere. 



2. OSSE AND RXre SPECTRA 

Table 1 gives the log of the 9 OSSE observations, together 
with results of power-law fits and basic data about the corre- 
sponding X-ray and radio states. The OSSE instrument ac- 
cumulated spectra in a sequence of 2-min. measurements of 
the source field alternated with 2-min., offset-pointed mea- 
surements of background. The background spectrum for each 
source field was derived bin-by-bin with a quadratic interpola- 
tion in time of the nearest background fields (see Johnson et 
al. 1993). Figure 1 shows the OSSE spectra (including stan- 
dard energy-dependent systematic errors), which were fitted up 
to energies at which the source signal was still detected. The 
uncertainty for a fitted parameter corresponds hereafter to 90% 
confidence (Ax^ = 2.71). We see that the source went through 
wide ranges of radio and X-ray fluxes and types of X-ray vari- 
ability during those observations. In spite of that variety, 8 out 
of 9 OSSE spectra are best-fitted by a power law with a photon 
index ofF ~ 3.0±0.1 and the flux varying within a factor of 2. 
The only exception is the OSSE spectrum corresponding to the 
highest X-ray flux measured by the ASM (1999 April 21-27), 
which is much harder, T ~ 2.3, and has a much lower flux. 

We then consider the OSSE spectra corresponding to the 
extreme X-ray fluxes measured by the RXTE/ASM, i.e., from 
1997 May 14-20 (VP 619) and 1999 April 21-27 (VP 813). 
We fit them together with spectra from the pointed RXTE ob- 
servations of 1997 May 15 and 1999 April 23 (the observation 
IDs are 20187-02-02-00,40403-01-07-00; 1% systematic eiTor 
is added to the PCA data with the responses of 2001 February). 
These PCA data correspond to the variability classes (Belloni 
et al. 2000) of x and 7, in which the variability is moderate and 
the source spends most of the time in two basic low (C) and 
high (B) X-ray flux state, respectively. 

We fit the data with the XSPEC (Arnaud 1996) model 
eqpair (Coppi 1999; G99), which calculates self-consistently 
microscopic processes in a hot plasma with electron accelera- 
tion at a power law rate with an index, Finj, in a background 
thermal plasma with a Thomson optical depth of ionization 
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electrons, Tj. The electron temperature, kT, is calculated from 
the balance of Compton and Coulomb energy exchange, as well 
as pair production (yielding the total optical depth of r > ri) 
is taken into account. The last two processes depend on the 
plasma compactness, £ = Ca^/ {TZnicC^), where £ is a power 
supplied to the hot plasma, TZ is its characteristic size, and ctt 
is the Thomson cross section. We then define a hard compact- 
ness, £\i, corresponding to the power supplied to the electrons, 
and a soft compactness, £s, corresponding to the power in soft 
seed photons irradiating the plasma (which are assumed to be 
emitted by a blackbody disk with the maximum temperature, 
fcTbb)- The compactnesses corresponding to the electron ac- 
celeration and to a direct heating (i.e., in addition to Coulomb 
energy exchange with non-thermal and Compton heating) 
of the thermal are denoted as £nth and ^tii' respectively, and 
£h = ^nth + £th- Details of the model are given in G99. 
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Fig. 1 .— (a) The 9 OSSE spectra of GRS 1915+105 deconvolved assuming 
power-law models. The color-coded spectra are labeled by the OSSE obser- 
vation numbers (Table 1) ordered according to the 50-60 keV flux, (b) The 
average spectrum of all observations with its power-law best fit. When fitted 
by an e-folded power law, the lower limit on the e-folding energy is 0.9 MeV. 

We also take into account Compton reflection with a solid an- 
gle of fl (Magdziarz & Zdziarski 1995) and an Fe Ka emission 
from an accretion disk assumed to extend down to lOGM/c^ 
(which results in a relativistic smearing). The equivalent width, 
Wko, with respect to the scattered spectrum only is tied to 
via WKa ^ 100(f7/27r) eV (George & Fabian 1991). The ele- 
mental abundances of Anders & Ebihara (1982), an absorbing 
column of > 1.75 x 10^^ cm'^ (Dickey & Lockman 1990; 
Vilhu et al. 2001), and an inclination of i = 70° are assumed. 

As discussed in G99, depends weakly on £s in a wide 
range of this parameter. An increase of is leads to increasing 
pair production, which then leads to an annihilation feature 
around 511 ke V. The presence of such a feature is compatible 
with the OSSE data (Fig. 1), but only very weakly constrained. 
G99 found that 4 = 10 provides a good fit to Cyg X-1 data. 
Here, we find that a good fit is provided with £s = 100, compat- 



ible with the high luminosity of GRS 1915-1-105. For example, 
for 1/2 of the Eddington luminosity, Le, and spherical geome- 
try, the size of the plasma corresponds then to ~ 100GA//c'^. 

This model provides very good description of our two broad- 
band spectra (as well as of other RXTE-OSSE spectra, S. V. 
Vadawale et al., in preparation). For VP 619, we assume a free 
relative normalization of the HEXTE and OSSE spectra with 
respect to those from the PCA. On the other hand, the HEXTE 
spectrum for VP 813 has relatively few counts at its highest en- 
ergies, and thus we use the actual OSSE normalization in that 
fit. Table 2 gives the fit results, and Figure 2 shows the spectra. 
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Fig. 2.— Fits to simultaneous PCA-HEXTE-OSSE spectra from VP 619 and 
813 with the hybrid Comptonization model. The dashed and solid curves show 
the models of the observed spectra and the unabsorbed spectra, respectively. 
The data are normalized to the PCA. 

Our model predicts the power law emission extending with 
no cutoff well above 1 MeV and a weak annihilation fea- 
ture (with the plasma allowed to be pair-dominated, i.e., with 
Tj 0, for VP 619). Those predictions can be tested by future 
soft 7-ray detectors more sensitive than the OSSE. We note that 
the COMPTEL has already detected a power law tail up to ^5- 
10 MeV in the soft state of Cyg X-1 (McConnell et al. 2000). 

Figure 3a shows the spectral components of the fit to the VP- 
619 spectrum. Compton reflection with i7 ~ 27r is detected at 
a very high significance (51 = gives x^/^^ = 214/131, cor- 
responding to the chance appearance of reflection of 6 x 10~^^ 
from the i^-test), and it is responsible for the convex curvature 
in the ^^10-100 keV spectrum. Figure 3a also shows that the 
scattered component has a spectral break at 100 keV but con- 
tinues as a power law (with addition of the broad annihilation 
feature) above it due to the domination of non-thermal scat- 
tering at those energies. Comptonization by the thermal elec- 
trons dominates at energies close to the blackbody component 
and thus PCA data of GRS 1915H-105 can be reasonably mod- 
eled up to 60 keV by thermal Comptonization of a disk black- 
body (Vilhu et al. 2001). When the OSSE data are included, 
the probability that non-thermal electrons are not present (i.e., 
i^jjjij — 0) is only 4 x 10^^". The best-fit thermal Compton 
model shown by the long dashes in Figure 3a strongly under- 
estimates the flux above 100 keV. The statistical significance of 
the presence of non-thermal electrons can be further increased 
by fitting the RXTE spectrum together with the average spec- 
trum from OSSE, which has virtually identical shape to that of 
VP 619, but much better statistics. Then, allowing for ^^th > 
reduces x^/^ from 228/139 to 105/137, which corresponds to 
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the chance probability of 4 x 10^^^. Thus, we strongly rule out 
the pure-thermal Comptonization model. 

During the VP 813, Compton reflection is statistically not 
required, as indeed expected at the large r ^ 4 of the scatter- 
ing medium covering the disk (which would completely smear 
out any disk reflection and fluorescence features). Thus, we set 

= is the fit. The presence of non-thermal electrons is now 
required at an extremely high significance (1 — lO^""^) due to 
the presence of the very distinct hard high-energy tail above the 
thermally-cut-off spectrum, see Figures 2, 3b. 
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Fig. 3. — (a) Components of the fit to the VP 619 data. All spectra are intrin- 
sic, i.e., corrected for absorption. The dotted, dot-dashed and dashed curves 
give the unscattered blackbody component, the scattered spectrum, and the 
component due to Compton reflection and Fe Kq fluorescence, respectively. 
The solid curve is the total spectrum. The thin long-dashed curve shows the 
best-flt thermal Comptonization model, which lies much below the data above 
100 keV. (b) The total model spectrum and the corresponding two components 
for the VP 813 data. The cutoff at > 10 MeV is due to pair absorption. 

On the other hand, Shrader & Titarchuk (1998) have fitted 
a model of bulk-motion Comptonization of blackbody photons 
to RXTE and BATSE data from a hard state of GRS 1915H-105 
similar to that of VP 619. We fit their model (bmc in XSPEC) at 
a free A'h to our broad-band spectra, and find it is completely 
unacceptable statistically, with 1^ = 851/133 and 549/92 
for VP 619 and 813, respectively. 

However, the specific feature of bulk-motion Comptoniza- 
tion is a high-energy cutoff at > 100 keV due to the effects 
of Compton recoil and gravitational redshift close to the black- 
hole horizon (e.g. Laurent & Titarchuk 1999). Such a cutoff 
is not included in the bmc model (and its inclusion would fur- 
ther worsen the fits above). Thus, we use Monte Carlo results 
of Laurent & Titarchuk (1999), which include the cutoff, to 
test whether the OSSE data (regardless of the data at lower 
energies) are compatible with its presence. We find that their 
theoretical spectrum for the accretion rate of M = 2L-e/c^ 
matches well the slope of the average OSSE spectrum at low 
energies (the histogram in Fig. 4). The Monte-Carlo spectrum 



can then be very well reproduced by a power law times a step 
function convolved with a Gaussian (model plabs ( step) in 
XSPEC, with the cutoff energy of 150 keV and the Gaussian 
width of 35 keV, the solid curve in Fig. 4). We see that the 
OSSE average spectrum lies well above that model at > 100 
keV. Quantitatively, the bulk-motion Compton model yields 
X^/i^ = 745/48. In comparison, the power-law and eqpair 
models yield /v = 29/48 and /v = 31/44, respectively. 
Thus, the bulk-Compton model is completely ruled out. Further 
problems with that model are discussed in Zdziarski (2000). 




Fig. 4. — Comparison of the average OSSE spectrum fitted by a power law 
(crosses) with predictions of the bulk-motion Comptonization model (Laurent 
& Titarchuk 1999; red histogram). The blue solid curve is an analytical ap- 
proximation of that model, see text. The theoretical spectrum predicts a correct 
low-energy slope, but it fails to reproduce the data at > 100 keV. 

Also, the commonly used phenomenological models of disk 
blackbody and either a power law or an e-folded power law give 
very bad fits to our data. The latter yields x^ jv — 263/132, 
405/91 for VP 619, 813, respectively. In fact, even the thermal 
part of the VP-813 spectrum is very poorly described by a disk 
blackbody, with -^jv = 3450/36 for a fit to the 3.5-20 keV 
PCA data, whereas the same data are well modeled by thermal 
Comptonization, x^ jv — 22/34 (with neither model includ- 
ing a high-energy tail). Thus, we find physical models of the 
spectra of GRS 1915H-105 in terms of thermal and non-thermal 
Comptonization and (in some cases) Compton reflection to be 
vastly superior to any other model proposed so far. 

3. DISCUSSION AND CONCLUSIONS 

We have found that broad-band spectra of GRS 1915H-105 in 
its two main spectral states (B, C) are very well fitted by Comp- 
tonization of disk blackbody photons in a (hybrid) plasma with 
both electron heating and acceleration. The presence of strong 
reflection indicates the plasma is located in coronal regions 
(possibly magnetic flares) above the disk. This physical model 
is the same as that fitted to the soft state of Cyg X-1 by G99. 
Differences in the variability properties of Cyg X-1 and GRS 
1915H-105 are likely to be due to the disk being stable in the for- 
mer case and unstable in the latter, most likely due to its much 
higher AI. The corona is unstable in both cases. The Comp- 
tonizing medium in the high-L state (B) is Thomson-thick, and 
it can represent the surface layer of an overheated disk accret- 
ing at a super-Eddington rate (Beloborodov 1998). An issue 
to be addressed by future research is the origin of the harden- 
ing of the tail in this state as compared to other ones (Fig. 1). 
Our model predicts broad annihilation features, although their 
strength depends on the unknown size of the plasma. 
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GRS 1915+105 was in the power-law 7-ray state in the clas- 
sification of Grove et al. (1998) during the 9 OSSE obser- 
vations. This state usually corresponds to the high/soft X- 
ray state. Indeed, X-ray spectra observed so far from GRS 
1915+105 (VUhu et al. 2001) are substantially softer than those 
with r ~ 1.7 and a sharp thermal cutoff at > 100 keV, char- 
acteristic to the hard state of other black-hole binaries. 
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TABLE 1 

OSSE OBSiiR\,\riONS Ol- GRS 1915+105 



OSSE 


Start 


End 


Exposure*" 


r 




OSSE 


X-ray 


X-ray flux" 


2.2 GHz flux 


ypa 


date 


date 


[10^ s] 






flux'= 


obs.'l 


[mCrab] 


[mJy] 


507 


1995 Nov 28 


Dec 7 


2.26 


3.081 


o.os 

0.07 


I8.8l";4 








601 


1996 Oct 15 


Oct 29 


7.15 


3.081 


0.05 
0.05 


-'-■-'--0.2 




840±330f 




619 


1997 May 14 


May 20 


3.27 


3.061 


0.06 
0.06 


-■-^■'-0.4 


axxotx 


270±120 


33±12f 


720 


1998 May 5 


May 15 


4.20 


3.OOI 


0.09 
0.09 




xxxxx 


570±40 


91±22 


811 


1999 Apr 6 


Apr 13 


2.92 


3.I0I 


0.07 
0.07 


14.9l'>-t 


app 


830±350 


11±6 


812 


1999 Apr 14 


Apr 20 


2.85 


2.891 


0.10 
0.10 


5+0-4 




920±270 


9±4 


813 


1999 Apr 21 


Apr 27 


3.23 


2.33I 


0.28 
0.27 


3 0+0-4 


1 


1190±240 


10±4 


917 


2000 Apr 18 


Apr 25 


1.56 


3.161 


O.OS 
0.07 


-'-^■-'--0.4 


xxxx<y-o- 


430±100 


19it7 


919 


2000 May 9 


May 26 


3.89 


2.961 


0.08 
0.08 


J^J^-O_0.3 


XPPPPPPPP 


700±90 




Sum 






29.60 


O.Ud_Q 04 


12 2+°-l 
'■■^■■^-0.2 









^ CGRO Viewing Period. 
Scaled to a single OSSE detector 
For 50-300 keV in units of lO^^o erg cm 5-1. 

Pointed RXTE observations with Greek letters giving their variability in the classification of Belloni et al. (2000). 
The RXTE/ ASM 2-12 keV flux. 
' These errors give the standard deviation of the flux variability. 

TABLE 2 

Model parameters of the broad-band spectra 



OSSE 

VP 




rkeVl 


4/4 




Tinj 






[keVl 


f^/27r 








619 
813 


9+0.3 
6 1+°-^ 

°-^-0.4 


1 07+0.05 
^•■^'-0.05 
, ,c+0.04 
1.00_Q 02 


n 4S+0-03 
^■4«-0.03 

p 90+0.02 


12+°°2 


9 9+0.2 


4 42+0.14 


0.39 
4.43 


66 
3.6 


9fi+0-50 
U.»D_Q 26 

(fixed) 


3.5 
8.9 


6.5 
17 


91/130 
60/90 



^ The total optical depth and the electron temperature calculated from energy and pair balance for the best-fit model (i.e., not free parameters). 

*" The unabsorbed bolometric flux of the model in units of 10^** ergcm"'^ s~^. 
The model luminosity assuming isotropy and a distance of 12.5 kpc in units of IO'"** ergs^^. 
The values oty^ jv < 1 result here from the assumed systematic errors, and are not due to overmodeling. 



